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PACS 12 . 20 . -m - Quantum electrodynamics 

PACS 44 . 40 . +a - Thermal radiation 

PACS 42 . 25 . Fx - Difltraction and scattering 

Abstract. - We discuss non-equilibrium extensions of the Casimir force (due to electromagnetic 
fluctuations), where the objects as well as the environment are held at different temperatures. 
While the formalism we develop is quite general, we focus on a sphere in front of a plate, as 
well as two spheres, when the radius is small compared to separation and thermal wavelengths. 
In this limit the forces can be expressed analytically in terms of the lowest order multipoles, 
and corroborated with results obtained by diluting parallel plates of vanishing thickness. Non- 
equilibrium forces are generally stronger than their equilibrium counterpart, and may oscillate 
with separation (at a scale set by material resonances). For both geometries we obtain stable 
points of zero net force, while two spheres may have equal forces in magnitude and direction 
resulting in a self-propelling state. 



The original quantum Casimir effect [T] is due to zero 
point fluctuations of the electromagnetic (EM) field in 
the vacuum between perfectly reflecting objects. Subse- 
, quently Lifshitz [2J treated the more realistic case of di- 
electric media at finite temperature by considering fiuc- 
tuating currents inside the objects, including both zero 
point and thermal fiuctuations. In general, the former 
dominate the force at small separations, while at separa- 
tions large compared to the thermal wavelength At, ther- 
mal effects prevail [21l3]. In situations out of equilibrium, 
the current fiuctuations in each body have to be treated 
separately at the corresponding temperature, e.g., using 
fluctuational electrodynamics introduced by Rytov over 
60 years ago [4]. Recently, out of equilibrium Casimir 
forces have been computed in a number of cases including 
parallel plates [5], modulated plates [5], as well as a plate 
and an atom in different setups There also exists 

a large literature on forces between atoms or molecules in 
non-equilibrium [lOHlSj . Formalisms for treating multi- 
ple objects at different temperatures have been recently 
presented [14,15^. In particular, for compact objects, ra- 
diation from the environment contributes to the force and 
has to be incorporated. 

Here, we treat (analytically as well as numerically) the 



cases of two spheres and a sphere in front of a plate. 
Keeping the description as simple and concise as possi- 
ble, we focus on the regime where the spheres are small 
compared to the separation (non-equilibrium effects are in 
most cases negligible at small separations), as well as ther- 
mal wavelengths. These restrictions allow the use of a one 
reflection approximation, as well as limiting to the spheres' 
(frequency-dependent) dipole response, respectively. We 
flnd a variety of interesting effects: The forces can be re- 
pulsive, oscillate or admit stable (zero force) points. At 
large separations, non-equilibrium forces decay as 1/d^ for 
two spheres and become independent of distance for sphere 
and plate. We also flnd points in which a pair of spheres 
experiences forces of equal magnitude in the same direc- 
tion. In the absence of other forces, this leads to a coop- 
erative motion of two identical spheres at constant sepa- 
ration, i.e. a self-propelled state. There are similarities to 
studies of atoms in non-equilibrium situations which we 
shall briefly comment upon. 

As presented in Ref. [15], our formalism treats N ob- 
jects (labeled as j = 1 . . . A^) in vacuum, held at constant 
temperatures {Tj\, and embedded in an environment at 
temperature Te„„. The conceptual starting point is the 
EM fleld radiated by isolated objects, each at its respective 
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temperature, which is then scattered by all objects while 
the environment is treated as an additional embedding 
"object." Physical quantities are then computed from the 
correlation function C'"'^'' of the electric field E at fre- 
quency Lo and points r and r' (both outside all objects) 



C"^'(Te„„,{r,}) = (E(c.;r)®E*(c.;r')) 

C'^iTenv) + ['^T^Tj) - Cf (Te„,)] . (1) 

3 



Equation ([TJ highlights the contribution of the dif- 
ferent temperatures to the non-equilibrium correlation: 
C'"^{Tenv) is the equilibrium correlation, i.e., with all tem- 
peratures held at T^nv (and including zero point fluctua- 
tions). It leads to the equilibrium Casimir force at tem- 
perature Tenv and is regarded as known. The difference of 
C'^'^'iT^nvATj}) from C"=«(re„„) is due to the deviations 
of the object temperatures Tj from T^nv Although deal- 
ing with -I- 1 sources, we have thus only to evaluate the 
N terms {Cj'^(r)}, the field correlations sourced by object 
j and scattered by all objects. In Ref. [15], we showed that 
Cj'^{T) can be derived by first considering the radiation 
of the object in isolation. 



vj Im Sj G* , 



(2) 



where axiuj) 



■{ex.p[hu!/kBT] — 1) ^, and Gj is 



the Green's function of the object. The QED origin of the 
force is manifested by the speed of light c and Planck's 
constant h. Cj{Tj) is found by integration over the envi- 
ronment sources (TSldS], subsequent scatterings lead to 



Cf(T,) ^ 0,Q(T,)0], 



with 
1 



(3) 



The multiple scattering operator Oj is expressed in terms 
of the composite T-operator Tj describing scattering by 
the other objects (not j), as well as the free Green's func- 
tion Go. For two objects, Tj = T2 is the operator of the 
second object. The force F acting on one of the objects 
(say object k) in this non-equilibrium situation is given 
by the integration of the Maxwell stress tensor cr over a 
surface Sk enclosing only this object, projected onto the 
surface outward normal n^. 



F'' = Rc i a- rifc dA . 



(4) 



The stress tensor is related to the field correlations, since 

'Jabir) = (e^EI + B^Bl - \ {\E\^ + |S|2) <5,,^ , 

where a, 6 = 1,2,3. Note that the sum of forces on all 
objects does not necessarily vanish, (i.e., there can be a 
net force on the system), and we must consider the force 



acting on each object separately. From Eq. ([T]), F*^ has 
the following contributions 



F'=(T_,{r,}) = F'=-^^(r_) ^-Y^\e){T.^ - f^^(t_)] 



(5) 



Here, F*''''^'(Te„„) is the force in equilibrium, and Fj^(Tj) 
is the force acting on object k due to the sources in object 
J at temperature Tj (obtained from the stress tensor in 
Eq. (HI) for the field Cf (T,)). 

Let us first consider two spheres of radii Rj (j = 1,2) 
with complex dielectric and magnetic permeabilities Sj 
and /ij, at center-to-center distance d and temperatures 
Tj, embedded in an environment at temperature Te„^. 
We derive the total force F^ acting on sphere 2; F^ is 
then found by interchanging indices 1 and 2. In Eq. ([5]), 
F^ has three contributions: The equilibrium force for the 
two spheres evaluated at the temperature of the environ- 
ment, a contribution due to the deviation of Ti from Te„t, 
(Ff ) and a contribution due to the deviation of T2 from 
Tenv (F2). The force F^ follows from the heat radia- 
tion of sphere 1, which can be written in terms of its T- 
operator [TSlfTT]. For the case d » Rj considered here, a 
one reflection approximation for the operator in Eq. ([3]), 
Oi ~ (1 — G0T2), is asymptotically exact. It amounts to 
a one time scattering of the field radiated by sphere 1 at 
sphere 2, and subsequently performing the integration in 
Eq. ([4]) over a surface enclosing sphere 2. This integra- 
tion in terms of spherical waves has been discussed, e.g. 
in Ref. [TB]. The force F2 is calculated similarly, only 
here we consider the heat radiation of sphere 2, which is 
once scattered by sphere 1 , and the surface of integration 
closed around sphere 2. Consistent with symmetries, the 
force F^ in Eq.® is parallel to the axis connecting the 
spheres; we shall denote this component by F'^ and adapt 
the notation where a positive sign corresponds to attrac- 
tion. The resulting force [iS', contains the T-operators as 
well as translation matrices for spherical waves, organized 
in a series of multipoles of orders I. To terminate the se- 
ries at the dipole order (/ = 1), for the equilibrium force 
between spheres [20ll2T] , it is sufficient to require d^ Rj, 
while in the non-equilibrium case, we have to additionally 



require At ^ Rj (At 



7.6/xm at room temper- 



ature). This ensures R* = Rju/c <^ 1 for all relevant 
frequencies, and we restrict to terms linear in the two T- 
operators Tj^ = 7J^=i(w) for polarization P = N,M and 
/ = 1. Then F^{T) is 

lim = 



> du 



Re[r2^']+Im[r2^'] 



- 1 pp, 

ISc^ Slc'^ 



Re[r/ 



9c^ 



9c3 



u^d^ uj^d^ uj'^d'^ 



)ppi 



(6) 



^We note that Ref. |14l performs a different decomposition, in- 
volving the equihbrium force at the temperatures of the objects 
rather than at T^nv as in Eq. JS]!. 
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In Eq. © (and Eqs. ©, ^ and (dH) below) we omit 
terms quadratic in Tj^ for brewity, a simplification which 
is justified for the cases considered beloA^0. 

For large separations, decays as and is repulsive. 
This originates from momentum transfer to the second 
sphere via absorption or scattering of photons. The other 
terms in Eq. ([6]), with higher powers in l/d, are (in most 
cases) attractive. Similarly, the seZ/ force F|(r), reads 



hm 

{d,AT}»-Rj 



h 

CTT 



duj- 



UJ 



9c' 



e''BT - 1 



.^Re[r/ 

p 



+ 1- 



/162c6 .81c^ 



(7) 



originating from radiation of sphere 2. Here P — M if 
P — N and vice versa. In contrast to F^, this term can 
oscillate as function of d at a scale set by material reso- 
nances. These oscillations originate from interference of 
two coherent traveling waves from sphere 2: i) a wave go- 
ing to sphere 1, being reflected back past sphere 2, and ii) a 
wave emitted in the reverse direction. Depending on wd/c, 
one has constructive or destructive interference. As these 
waves interfere in the exterior region, we expect the oscil- 
lations to become weaker as R2 becomes large compared 
to the penetration (skin) depths. For a sharp resonance 
of £2(w) at cjo in Eq. Q, the oscillations as function of 
distance have wavelength ttc/wq- 

For small insulating spheres of radius Rj (with /ij = 1), 
we employ the following expansions of the T-operator 



n-N -2^^ , , 



o 



(i?f) 



(i?f), (8) 



in terms of the complex frequency dependent dipole po- 
lariz ability. 



(9) 



Higher multipoles T^i for / > 2 are of order i?*^, and 
Eqs. dH) and can Idc simplified by use of Eq. The 
range of applicability of this approximation depends on 
material properties. An expansion of in both R* and 
y/eR* shows that the condition \y/e\Rj ^ At (in the rel- 
evant frequency range) is sufficient for many materials, 
including the ones studied below. For \ej\ >• 1 (conduc- 
tors), the expansion is generally not applicable (e.g. Tj^^ 
is then of order Rf [21]). With Eqs. JH]) and jH]), one sees 
that Fj is only nonzero if Im Sj 7^ (or for magnetic ma- 
terials Im/ij ^ 0), as only lossy spheres emit heat. This 
holds for any Rj. 



quency range 

|7^^P and vanishes for |e[ 



in the relevant fre- 
E.g. the emissivity of a sphere contains — ^3,^!■-^P^ 



^It requires {| Im[7^:P]|, | Re[7]' 



» 17-^1 = 



00 or Im[e] — > where Re [7^^] 
-|7^^|^, captured only by inclusion of the quadratic terms |19l . 
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Fig. 1: Total force on sphere 2 in a system of two Si02 spheres 
at separation d in a cold (0 K) environment. Dashed lines in- 
dicate repulsion. The crossing of solid red and dashed green 
curves represent a point where the forces are equal in mag- 
nitude and direction, see main text. Points of change from 
repulsive to attractive with increasing d are stable equilibria. 



The leading low temperature behavior of the force for 
insulators can be derived by requiring At S> Aq, where 
Ao is the wavelength of the lowest resonance of the mate- 
rial. The dielectric functions and polarizabilities are then 
expanded as [22] 



aj{uj) = ao,j + ioiio,' 



0(^2), 



(10) 

(11) 



with eo.j", Ainj, aoj and a^oj = 3i?^/(£o,j + 2)^ real. For 
At S> Aq , the interaction term is then given in closed form, 



lim Ff 



{ci,AT}>i?3 

+ 010,2 



5A7 



327r5AT Stt^A^ 



2\d 



5d3 



ISttA^ 



(12) 



The self force F| does not oscillate to lowest order in tem- 
perature and takes a lengthy form [19| . In the limit where 
d is the largest scale, we have 



lim 

d»AT>{fl3,Ao} 



Trd^ 



(13) 



While in this range of d the force F2 is independent of tem- 
perature, it vanishes as T — > since with At the largest 
scale (At >• {d, Rj,Xo}), one has 



hm F^ = 



ff)2 -^»",2QiO, 20^0,1; 



(14) 



which is identical to F^ in this limit, with indices 1 and 2 
interchanged. 

We evaluate the total force in dipole approximation nu- 
merically for i?i = i?2 = 1 /ini using Eq. ([5]) and the 
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Fig. 2: Total force on sphere 2 in a system of two Si02 spheres 
at separation d in a warm (3007^) environment. Dashed Unes 
indicate repulsion. The thin dotted line is the red line divided 
by 2, see main text. 
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Fig. 3: Total forces on Si02 spheres with shifted optic res- 
onances. Dashed lines indicate repulsion. Unstable (U) and 
stable (S) SPP and zero force points are visible. 



equilibrium Casimir Polder force (Eq. (94) in Ref. |23j) 
which in the relevant limits reads 



1- T^lea 161 /iC 

hm F ■ = -— -7^ao,iao,2, 

1- T^2 en l^^C 

d>{i?.j,Ar} a' At 



ao,i"o,2- 



(15a) 
(15b) 



Figure [T] shows the forces on Si02-spheres (we used op- 
tical data with Eq ~ 3.7) in a cold (0 K) environment. 
We evaluated Eqs. © and together with dH), 
Within these simplifications, the forces are proportional 
to R\R\ {R~ I ^m is roughly the upper bound of valid- 
ity of this approximation for Si02 at room temperature, 
where for the total heat emitted by an isolated sphere, the 
asymptote a differs by 12% from the exact result [T5]). 
The force starts to deviate strongly from its equilibrium 
value around d « At/2. Sphere 2 is repelled at large d if 
Ti = 300 JsT due to the radiation pressure. If additionally 
T2 = 300K, the oscillating force F| is visible and it dom- 
inates the total force for large c? if Ti = OK; the net force 
now has many zero crossings, where every second one is 
a stable equilibrium point. As discussed above, we expect 
the wavelength of the oscillations to be roughly 4.75 /im 
due to the resonance of Si02 at wavelength 9.5 /im. Addi- 
tional modulations are due to interferences with a second 
resonance of Si02 at 22 /im. The figure also provides com- 
plete information about the force on sphere 1: e.g., in case 
Ti = and T2 = 300 K, the red curve shows the force act- 
ing on sphere 2, while the green curve shows the force 
on sphere 1. At the crossing of the solid red and dashed 
green curves the two spheres feel equal forces in the same 
direction. This corresponds to what we define as a self- 
propelled pair (SPP), where the spheres experience equal 
acceleration in the same direction and hence remain at a 
fixed separation. Note, however, that this is an unstable 
arrangement in which any small perturbation leads to the 
spheres moving apart. 

Figure [5] shows the situation for a warm {300K) en- 



vironment. Here, the force has repulsive parts only if 
Ti = 300K and T2 = OK where it shows multiple stable 
equilibrium points. For all other cases, the force is purely 
attractive, decaying as l/d^ if Ti = 0. For Ti — and 
T2 — 300K, one has stable and unstable SPP's, e.g., where 
the black dotted curve crosses the dashed green curve in 
Fig. [21 i.e., for R2 — i?i2^/'^ assuming solid spheres with 
mass ^ R^. 

The l/d^ contribution to in Eq. ^ (with Eq.®) 
is proportional to the product of the imaginary parts of 
the polarizabilities. These are peaked at the resonances of 
the material and this term can be suppressed by reducing 
the overlap of resonances. Figure [3] shows the forces where 
the dielectric function of sphere 2 is replaced by e2{(^) = 
esi02 (l-l^w), which, in principle, can be achieved by using 
different isotopes. Now, the forces are asymmetric even for 
Ti = T2, and due to the suppression of F^, we have stable 
as weU as unstable SPP's for e.g. Ti = = 300K, and 
Tenv = OK, in contrast to Fig. [TJ 

For a sphere (i?, Ss, /i^, T,) in front of a plate (Sp, 
IJp, Tp) at center to surface separation d, Eq. (O gives 
distinct non-equilibrium forces acting on the plate, or on 
the sphere. While both can be derived with equal effort, 
we restrict to the force acting on the sphere, = — • 
Hp (with outward normal Hp of the plate), separated into 
Fp and Fj*. Scattering from the plate is governed by the 
Fresnel reflection coefficients r^ for P = M, N, given by 



with r^ obtained from r*^ by interchanging and e. In 
the one reflection approximation, the force F^ is derived 
by a one-time-scattering of the radiation of the plate at 
the sphere. The subsequent integration in Eq. ^ is done 
in plane waves basis, over two planes, enclosing the sphere 
and separating it from the plate. For the contribution 
F/, the sphere radiation is scattered at the plate with 
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identical surface of integration. As before, this procedure 
(valid for R) 19 involves all multipoles of the sphere. 
Only for i? ^ At, we can further restrict the T operators 
to I = 1. The interaction term Fp{T) has two distinct 
contributions, Fp = F^ ^^ + F^ ' 
evanescent waves emitted by the plate 



p,pr • -p,eDJ from propagating and 



lim F' 

{d,AT}»i?. ^ 



2c7r 



duj- 



ifpr + fev), (16) 



e ''B ' 



where the functions 

2 /•'^/c 



k I dk 







E 

p.p' 



fev — 2 I 



k I dk 



±aK±e 



(1 - ir 



Re[r 



p'l 



(17) 



2/^2 



oj/c 



p 



Im 



Im[r^ 



(18) 



explicitly contain the radiation of the plate l4l. The force 
Fp is d independent as it arises from absorption or scat- 
tering of far field photons by the sphere, while the near 
field contribution Fp depends on d. The self term (T) 



lim F' 

{d,AT}>fl 



TT 



P 



duj- 







1) 







Re 



A:2 c2 



- 1 



k±dkj 



(19) 



contains both evanescent and propagating contributions 
but no separation independent term. Instead, F^ behaves 
similarly as F| in Eq. ([7]), oscillating as a function of d, 
falling off at large separations as 1/d. 

For a dielectric sphere and plate, we next employ 
Eqs. pO)) and (fTTj) . to obtain the leading behavior at low 
temperatures (Ay ^ {Ao,i?}, but not necessarily d). The 
d independent part now becomes. 



lim F' 



"63" 4^^"^'^ 



0)A,„,sa,o. (20) 



- p.ev 



can be analyzed in the following two limits, corre- 



sponding to expansions of the function fev{i^,d), 



lim 



F' 



TT he 



6 A|d3 



Re 



1 + £o,i 



1 



ao. (21) 



<i»Ar>{fl.Ao} 

In the opposite limit, with At 3> {d, -Rj, Aq}, we have 

1 



2 A2,d4 (i+^„_p)2-«o- (22) 

Equation (HT)) is similar to Eq. (12) in Ref. [5]. As was the 
case for F| , in leading order in temperature the self part 
F/ does not oscillate. For d » At ^ {i?, Aq}, we have 
FJ* c>c l/d®, the counterpart of Eq. ([T^ . with a lengthy 
prefactor. For At ^ {d, Rj, Aq} we have 



,• TPS eo,p - 1 



4 A2 d4 eo.p + 1 ' 



(23) 




d [urn] 

Fig. 4: Total force on a Si02 sphere of i? = \\xra. in front of a 
Si02 plate in a cold environment. Dashed lines indicate repul- 
sion. Every second zero of the red curve is a stable equilibrium 
point. 



when interchanging real and 
:. The equilibrium force can 



which is identical to Eq. (122 



imaginary parts for 
be found in Ref. [21 



and 

For d 3> i? one has 



lim F^'^« = ^5^^^ao$(£op), 

lim F ' ' = —7- ao, 

d->{R,\T} Ad^XT £o,p + 1 



(24a) 
(24b) 



where $(eo.p) is e.g. given in Ref. [24]. Figure |4] shows 
numerical results for the force on a sphere in front a plate 
(both made of Si02) for R = 1/xm in a cold (0 K) environ- 
ment. Again, we use the simplification of Eq. ([5]) and the 
resulting force is proportional to R^ (also here, R — 1/xm is 
roughly the upper bound of validity of this simplification). 
If the plate is warm, the distance independent repulsion 
is visible. If only Tg is different from Te„^, the force F^ 
dominates at large d, leading to multiple stable points. 

Figure [S] shows the curves for a warm (SOOJC) environ- 
ment. Here, the d independent force (for Tp = 0) is attrac- 
tive. Again, if only is different from Tenu, we observe 
many changes of the sign of the force. Exploring the ef- 
fects of shifting resonances, we found that in contrast to 
the case of two spheres, here shifting suppresses the self 
term more strongly than the interaction term. For the 
special case of a (resonance-shifted) SiC spher^ in front 
of a Si02 plate, see Fig. [51 the temperature of the sphere 
is almost irrelevant for the force. This is beneficial to ex- 
perimental setups, as it is presumably harder to maintain 
Tg at a constant value, compared to keeping Tp and T^nv 
constant. Additionally, in Fig. [51 the special choice of pa- 
rameters leads to a stable equilibrium point, which, again, 
is almost independent of Ts . 

The presented formulae for the forces, i.e., Eqs. (O, ©, 



T6l) and (HH) after substitution of F 



N 



:2ui- 



a 



j{uj) and 



^We found that the described effect appears most pronounced by 
adjusting the resonance of SiC (see Ref. [25) for s{u})) by insertion 
of a factor of 0.75 in frequency. 
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Fig. 5: Total force on a Si02-sphere of i? = l/im in front of 
a Si02 plate in a warm {300K) environment. Dashed lines 
represent repulsion. 

the magnetic dipole polarizability, can be derived inde- 
pendently from diluting two plates of vanishing thickness, 
confirming the correctness of our formalism for compact 
objects. This calculation will be presented elsewhere. 

Our results constitute a macroscopic generalization of 
non-equilibrium interactions between thermal gases of 
atoms and interactions between atoms in excited states 
0. However, we emphasize that the forces on two macro- 
scopic objects are not equal and and opposite, an effect 
which cannot be found from the interaction potential as 
used in studies of two atoms [TUHS]- We hope that our 
results may eventually shed new light on the debated non- 
equilibrium interactions of atoms. 

While, for simplicity, we discussed the forces for small 
radii and moderate temperatures, our formalism is more 
generally applicable for any values of R, d and T. Future 
work will consider the cases of larger spheres where non- 
equilibrium effects may be stronger. 
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